Experimental detection of the diamino-pentazolium cation and theoretical exploration of derived high energy materials

In this work, we realized the detection of diamino-pentazolium cation (DAPZ+) in the reaction solution experimentally and proved it to be meta-diamino-pentazole based on the transition state theory. Quantum chemical methods were used to predict its spectral properties, charge distribution, stability and aromaticity. Considering that DAPZ+ has excellent detonation properties, it was further explored by assembling it with N5−, N3− and C(NO2)3− anions, respectively. The results show a strong interaction between DAPZ+ and the three anions, which will have a positive effect on its stability. Thanks to the high enthalpy of formation and density, the calculated detonation properties of the three systems are exciting, especially [DAPZ+][N5−] (D: 10,016 m·s−1; P: 37.94 GPa), whose actual detonation velocity may very likely exceed CL-20 (D: 9773 m·s−1). There is no doubt that this work will become the precursor for the theoretical exploration of new polynitrogen ionic compounds.

NMR data.Unfortunately, the half-life of the aminopentazole derivative in d 7 -DMF solution was estimated to be around only 11 min at 21 °C.Very recently, Lu and Xu et al. attempted to use Oxone to oxidize cyclo-N 5 − to obtain covalent derivatives N 5 OH (4) or N 5 O −25 .However, due to a very high activation energy barrier involved in the reaction, they were unable to obtain the expected pentazolate N-oxide.Overall, despite the tremendous efforts of researchers from around the world, the syntheses of covalent pentazole derivatives remain a significant challenge due to their poor stabilities.
In our continued efforts to explore pentazole derivatives, we tend to use cyclo-N 5 -ionic derivatives as precursors for the construction of covalent bonds on the pentagonal ring.Fortunately, it is proud that in this work, based on our previous research, the diamino-pentazolium cation (DAPZ + ) was successfully prepared and characterized, and the related contents were supplemented by the theoretical calculation method of quantum chemistry.Considering its high energy, the assembled ionic compound system was analyzed in detail by Atoms in Molecules (AIM), Independent Gradient Model based on Hirshfeld Partition (IGMH), Electrostatic Potential (ESP) and Symmetry-Adapted Perturbation Theory (SAPT), and the detonation performance was also taken into account. 6in hand, we focused on the corresponding N-amination reaction using O-p-toluenesulfonylhydroxylamine (THA) and HOSA.However, the N-amino product could be prepared only by HOSA instead of THA.

Synthesis and reaction mechanism of m-DAPZ
The amination of cyclo-N 5 − was carried out at 30, 40, and 50 °C, respectively, and an interesting phenomenon was found.The amination of cyclo-N 5 − was not successful at 30 °C, and the mass spectra showed that the amination reagent and the cyclo-N 5 − in the reaction mixture.At 40 °C, obvious peaks of a diamino-substituted pentazole species were observed at m/z = 100.945/102.128and mass-selected for subsequent MS/MS studies.No amination product was generated in the reaction system at 50 °C, and only the signal peaks of the amination reagents were present in the mass spectra, which may be related to the decomposition of cyclo-N 5 -precedes the reaction.Except for the signal of 100.945,HSO 4 − , and SO 3 2− were confirmed by secondary mass spectrometry to be by-products of amination reagents HOSA (Fig. S1).Further confirmation was obtained by MS/MS analysis of the m/z = 100.945peak in negative ion mode (Fig. 2a).Using a low collision voltage of 30 V, the m/z 100.945 ion was not shattered (Fig. 2b).As the collision voltage increased to 120 V, the m/z 100.945 ion underwent stepwise NH 2 loss, giving rise to an intense peak with m/z value of 84.950 (Fig. 2c).Then the voltage continued to increase to 200 V, the peak of m/z = 100.945almost disappeared and underwent stepwise second NH 2 loss rise to an intense peak with m/z value of 68.955 (Fig. 2d).According to the mass spectrum in positive ion mode, the peak of m/z = 102.128was detected, the m/z = 86.097and 70.066 were obtained by MS/MS analysis using a low collision voltage of 30 V which underwent twice NH 2 (Fig. 2e, f).
The pentazolate species replaced by two amino groups (DAPZ) cannot be identified from the mass spectrum as positive ion, negative ion or neutral molecule (Fig. S2).Therefore, different structures of DAPZ, DAPZ -, and DAPZ + were optimized by Gaussian09 at the M06-2X/6-31+G(d,p) level.However, it was impossible to calculate neutral DAPZ molecules (53 electrons) with multiplicity 1.For the negative ion state, the optimized configurations of the five DAPZ -ions obtained are not reasonable (Fig. S3).Finally, only the positive ion state is suitable.The five optimized DAPZ + structures having no imaginary frequencies and their total energies, enthalpies, free energies are shown in Fig. S4.Two of them are o-DAPZ + (o-1, o-2), and the other three are m-DAPZ + (m-1-m-3).We can find two pairs of chiral isomers with equal energies.According to the principle of minimum energy, the possible structure should be the bottom two m-DAPZ + (chiral isomers) in Fig. S4.
In order to further investigate the amination process and product structure in detail, the transition state of the whole process is calculated, as shown in Fig. 3, where IM refers to the intermediate substance and TS refers to the transition state found by different paths.Path I-III refers to pentazole amination, amino pentazole meta-amination and ortho-amination respectively, and the whole reaction path uses relative Gibbs free energy.It should be noted that the relative energies of the initial reactants of these three reaction paths are all set to 0, but their absolute Gibbs free energies are different, so as to compare the reaction barriers of different paths.The intrinsic reaction coordinate (IRC) calculations were performed on the details of the amination process, redisassembled and integrated the process in which the transition state reaches the reactant and the product from two directions, obtaining a complete reaction process from reactants to transition state and then to products.Please refer to the MP4 file in the supporting information.There are four videos provided, which correspond to the IRC processes of four reaction transition states, and more calculation details are in the following sections.
In the process of amination of cyclo-N 5 − , the potential barrier of TS1 is 27.68 kcal•mol −1 , and the Gibbs free energy of the intermediate product IM1, amino-pentazole, is lower than that of cyclo-N 5 − , which means that the amino-pentazole system is relatively more stable than that of cyclo-N 5 − .IM1, was further aminated, and the initial states of ortho and meta-amination were IM2 and IM3 respectively.Two different transition state paths show completely different information.From the distribution of transition states, the transition state of orthoamination TS2 is above the transition state of meta-amination TS3.It means that the barrier of ortho-amination will be higher and the reaction will be more difficult, which is confirmed by the comparison of Gibbs free energy.In addition, it is worth noting that in the process of ortho-amination, it can be seen that amino-pentazole cannot form m-DAPZ + in one step, but rather undergoes a hydrogen transfer process.However, the energy required for hydrogen transfer from the intermediate IM4 to o-DAPZ + is only 2.23 kcal•mol −1 , and the energy of the transition state of the reverse reaction is even less than 0.41 kcal•mol −1 , so o-DAPZ + will be easily converted into IM4.In other words, o-DAPZ + will be difficult to exist stably in the system.Another important phenomenon is that the final energy of m-DAPZ + will be lower than that of o-DAPZ + , which indicates that m-DAPZ + will be more stable in the system.After the above analysis, there is no doubt that most of the synthesized product is m-DAPZ + , but due to the technical limitations and the difficulty of synthesis, a very small amount of o-DAPZ + will probably exist in the whole system in the form of "impurities".
Positively charged m-2 (or m-3) is considered the most likely amination product, which has C 2 symmetry (Fig. S5).The bond lengths, bond angles, and torsion angles of m-2 are presented in Tables S1-S3.It consists of a five membered ring with a nearly planar structure and two meta-amino groups.The two amino groups are slightly deviated from the ring based on the torsion angles of 173.66° (N4-N3-N7-N11 and N8-N11-N7-N3).In the five-membered ring, the bond lengths of N1-N2, N1-N3 (or N2-N11) and N3-N7 (or N11-N7) are 1.305, 1.322 and 1.314 Å, respectively, which are shorter than that of the N-N (1.450 Å), implying they have properties of partial double bonds.The N-NH 2 bond (1.359 Å) is the longest N-N bond in the molecule, but it is also shorter than the conventional N-N bond.For the bond angles, m-2 is not regular polygon structures, and the bond angles are in the range of 100.33°-114.17°.
Atomic charge is one of the simplest and most intuitive forms of description of charge distribution in chemical systems and is important for the analysis of both chemical reactions and adsorption principles.Hirshfeld charge can be rectified using atomic dipole moment corrected Hirshfeld (ADCH charge), as the Hirshfeld charge levels are typically tiny and the dipole moment and electrostatic potential repeatability are quite low.In calculating the ADCH charge, the Hirshfeld charge and atomic dipole moment of each atom are first calculated, and then each atomic dipole moment is expanded into the corrected charge of the surrounding atoms, and sum them to obtain the ADCH charge.The charges of each atom in m-2 under the solvent model are shown in Fig. 4a.The charges on N3 and N11 of m-2 are more than those in N1, N2 and N7 because both N3 and N11 are connected to the electron-donating amino groups, which makes the charges on the amino group transfer to the adjacent nitrogen atoms on the ring.In addition, the nitrogen atom on the amino group transfers some of its charge to the hydrogen atoms, resulting in a reduction in its own charge, making the entire molecule more evenly charged and stable.Molecular electrostatic potential (ESP) of m-2 is analyzed by the Multiwfn program and shown in In order to provide experimentalists with spectral data for further instrumental identification of DAPZ + cation, we had calculated IR (Fig. 4c) and NMR (Fig. 4d) spectra of m-2.Obviously, the absorptions of m-2 are more complex than the cyclo-N 5 − , the new absorptions at 180, 706, 841, and 1654 cm −1 are attributed to the torsional vibrations of the N-H bonds.The strongest absorption at 3507 and 3628 cm −1 are attributed to the N-H stretching vibration of the amino groups.The N-N stretching vibration of the N-amino groups (N3-N4 and N8-N11) and the ring (N1-N3 and N2-N11) are located at 765 and 1540 cm −1 , respectively.The calculated 15 N NMR spectrum using nitromethane (Fig. S6) as TMS shows signals for the five-membered ring are at 32.008 (N1 and N2), 66.153 (N7), and 70.088 ppm (N11 and N3).In addition, the N4 and N8 of amino groups are observed at 286.233 ppm (Fig. 4d).
The HOMO-LUMO orbitals of m-2 and cyclo-N 5 − were also studied.The energy levels and the energy gaps are given in Fig. S7.The red phase means positive, while the blue is negative.The energy gap is an important value that can evaluate whether the electrons can be excited easily or not.In some ways, larger energy gaps mean higher chemical stability.As presented in Fig. S7, m-2 has much higher energy gap (12.85 eV) than the cyclo-N 5 − (11.82 eV), which indirectly indicates that m-2 is more chemical stable than cyclo-N 5 − .In addition, the delocalized electron clouds of m-2 and cyclo-N 5 − were visualized by local orbital localizer (LOL) analysis (Fig. 5a, b).In the π-electron density plane, the colors from dark blue to red show the low-high trend of π-electron density.Due to the influence of two amino groups, the aromaticity of m-2 is less than that of cyclo-N 5 − .Similar to cyclo-N 5 − , but smaller number of electrons gather on the five-membered ring, forming an almost closed electron ring in m-2 (Fig. 5c, d).Table S4 shows that the calculated bond dissociation energy (BDE) of m-2 is 342.82 kJ•mol −1 at 298.15 K.The large BDE indicates that once m-2 is synthesized, it will have a good stability, either spectroscopically or even in the condensed phase.
Finally, the thermal stability of DAPZ + must be mentioned.In fact, this is a common problem of all covalent pentazole derivatives, mainly because various types of chemical bonds formed between N atoms are easy to break and decompose under external stimulation, despite their high bond energy.Recently, a good approach was developed to estimate thermal stabilities of covalent and ionic pentazole derivatives using the longest N-N bond length in the ring 26 .The general conclusion is that the thermal decomposition temperature decreases with the increase of the longest N-N bond length.At the level of M06-2X/6-31+G(d,p), the bond length of the longest N-N in cyclo-N 5 − after optimization is 1.326 Å, while that in DAPZ + is 1.359 Å.It can be inferred that the thermal decomposition temperature of DAPZ + will be lower than that of cyclo-N 5 − .Further, it is not difficult to imagine that the synthesis and isolation of DAPZ + is very difficult.Fortunately, this instability can be alleviated by effectively enhancing the interaction between ions through hydrogen bonds, π bonds, and conjugated structures, thereby improving the stability of the whole system.Therefore, selecting appropriate anions is another part of this work.

AIM analysis
The theory of atoms in molecules (AIM) plays an important role in studying the interaction between atoms 27 .Topological analysis of electron density is a key part of AIM framework.The point in space with the first derivative of the vanishing electron density is called the critical point (CP).According to the sign of eigenvalues of electron density Hessian matrix, it can be divided into four categories, namely nuclear critical point (NCP), bond critical point (BCP), ring critical point (RCP) and cage critical point (CCP) 28,29 .Among them, BCP has attracted extensive attention because it can reflect the interaction between atoms.Figure 6 shows the electron density topology of three ionic compound systems, and the external diameters are also given 30 .Some key topological properties of these BCPs are listed in Table 1.
Generally speaking, the Laplace operator of electron density (∇ 2 ρ(r)) can reflect the type of interaction.When ∇ 2 ρ > 0, it indicates the existence of closed shell interaction (i.e.hydrogen bond and van der Waals interaction) 31 .Therefore, it can be clearly seen from the figure and table that there is interaction between some N and O atoms in the three systems.On the basis of Popelie's criterion for judging the existence of hydrogen bonds, Rozas further constructed a criterion for judging the strength of hydrogen bonds, that is, when ∇ 2 ρ > 0 and H(r) > 0, hydrogen bonds are weak; When ∇ 2 ρ > 0 and H(r) < 0, there is a moderate strength hydrogen bond.When ∇ 2 ρ < 0 and H(r) < 0, the hydrogen bond strength is high 32,33 .From the data in Table 1, it can be seen that H atoms on m-DAPZ + will form hydrogen bonds with anions in all three systems.The difference is that in the [m-DAPZ + ] [N 5 − ] system, the strength of hydrogen bonds formed is weak, for example, the electron density between H10… N16 is only 0.021612 a.u., its ∇ 2 ρ(r) is 0.066946, and H(r) is 0.000693, all of which are positive numbers.In the [m-DAPZ + ][C(NO 2 ) 3 − ] system, all hydrogen bonds with moderate strength are formed.These hydrogen bonds will play a very important role in stabilizing the system and enhancing the mutual attraction between anion and cation.But It must be admitted that the attractive interaction between m-DAPZ + and C(NO 2 ) 3 − will be stronger than that of N 5 − , and the system will be more stable.In addition, for the [m-DAPZ + ][N 3 − ] system, the first thing that can be determined is that weak hydrogen bonds must be formed between H10…N12.As for the interaction between H8…N14 atoms, although its atomic species conforms to the basic type of hydrogen bonding, its electronic density is 0.06311 a.u., which will not be within the scope of the existence standard of hydrogen bonding.It should be explained that the electron density between these two atoms is significantly higher than that between the atoms forming hydrogen bonds in the system, and there is no covalent bond between them, which means that not only hydrogen bonds may be formed between H8…N14, but also the strength of hydrogen Table 1.Electron density topological parameters at BCPs. www.nature.com/scientificreports/bonds will be higher.More importantly, although the electron density at the critical point of hydrogen bonding in the three systems is higher, and the interaction intensity of hydrogen bonding is stronger than that between other atoms (probably stronger van der Waals interaction), the interaction between these atoms also plays an important role in the stability of the complex.

IGMH and interaction energy analysis
Noncovalent interactions (NCI) analysis is a graphical method to study weak interaction, which can not only show the location of interaction directly, but also show the type and intensity of interaction at the same time 34,35 .Among them, by mapping the sign(λ 2 )ρ function with different colors, IGMH can not only clearly describe the types and properties of the interaction, but also flexibly customize the segments, and can only display the weak interaction between segments without interference from within the segments 36 .Therefore, in order to further explore the interaction between m-DAPZ + and other anions, the above three systems were analyzed by IGMH, as shown in Fig. 7. Different colors represent different types of interactions.Blue represents weak attractive interactions, such as hydrogen bonds and halogen bonds; Green represents van der Waals interaction; Red represents repulsive interactions, such as steric effects in cages and rings.First of all, in the system of [m-DAPZ + ][N 5 − ], it can be clearly seen that some N atoms on the N 5 -anion will form hydrogen bonds with H atoms on m-DAPZ + , which is consistent with the analysis results of AIM.In addition, in AIM analysis, there will be strong interaction between some N atoms on anion and cation, which can be fully explained in Fig. 7a, N 5 − and m-DAPZ + contain typical conjugated π bonds.When they form a composite system, they will inevitably show the form of conjugated π-π stacking, which explains that blue appears around every N atom in N 5 − in Fig. 7b.Secondly, other regions are mainly green, which means that the electron density is low, and there are extensive van der Waals interactions in this position.Among the three systems, the overall green region area of [m-DAPZ + ][N 5 − ] system is the largest.The difference is that in the [m-DAPZ + ][N 3 − ] ion system, in addition to the common hydrogen bonds, the N atom on N 3 − will still form visible interaction with the N atom on m-DAPZ + .Thanks to the special electronic properties of N 3 − , in Fig. 7c, the N atoms in the anion and cation will form LP (lone pair)-π interaction 37 , which is very common in the analysis of weak interaction of energetic materials.It is worth noting that these interactions play a key role in the mutual attraction of composite systems.An interesting phenomenon is that there are some red areas in both systems, which is obviously normal, because even in the cationic system, some atoms will be negatively charged and show polarity.Therefore, it is not difficult to draw a conclusion that anions will repel the more polar atoms on m-DAPZ + .In fact, the negatively charged N atom on the amino group in m-DAPZ + conforms to the above analysis and the position of the red area in Fig. 7 proves that this is the case.The IGMH of [m-DAPZ + ][C(NO 2 ) 3 − ] system is similar to that of the other two systems.m-DAPZ + will form strong hydrogen bonds with the active O atoms on it, and the van der Waals interaction between them will also stabilize the overall structure of [m-DAPZ + ][C(NO 2 ) 3 − ].The interaction energy of the three systems is decomposed and analyzed based on SAPT, as shown in Fig. 8. Generally speaking, among the three systems, electrostatic interaction plays a key role in stabilizing the composite system.Of course, although the contribution of dispersion and induction to the interaction between anions and cations is not as great as that of electrostatic interaction, it cannot be ignored.However, there is a negative correlation between the exchange mutual exclusion energy and the interaction between anion and cation.It is obvious from the figure that although the atom number of N 3 − is the least, its exchange mutual exclusion energy is the highest.This is mainly because the strong polarity of N 3 − has a great negative effect on the electron density between m-DAPZ + and N 3 − , leading to the weakening of the ion stabilization energy, which reflects the steric effects.Because of its conjugated structure, N 5 − has the lowest polarity, so the repulsion with m-DAPZ + is also the lowest.Nevertheless, some new conclusions can be drawn from calculating more accurate data.It has to be mentioned that although there are more hydrogen bonds or larger van der Waals interaction regions in the systems of [m-DAPZ + ][C(NO 2 ) 3 − ] and [m-DAPZ + ][N 5 − ], in fact, the stronger polarity of N 3 − will also bring some good aspects.From the overall interaction energy, the interaction energy of [m-DAPZ + ][N 3 − ] system is -530.59kJ•mol −1 , which is the smallest among the three systems.It is precisely because of the particularity of N 3 − electron system that electrostatic interaction is greater, and the mutual polarization and mutual transfer of charges between anion and cation are stronger (induced interaction), and their values are −584.11kJ•mol −1 and −107.21kJ•mol −1 respectively.It needs to be explained that this is the biggest as the exchange of mutual exclusion energy, and it is not in conflict.In contrast, the attraction interaction is more dominant.As for the systems of [m-DAPZ + ][N 5 − ] is stronger, which is also determined by the stacking of conjugated π-π.Therefore, it is not surprising that the [m-DAPZ + ][C(NO 2 ) 3 − ] ion system has the lowest interaction energy.

ESP analysis
Electrostatic potential (ESP) is widely used to study the intermolecular electrostatic interaction 38,39 and predict the reaction sites and molecular properties, so it is necessary to study the electrostatic potential of the system.Generally, the electrostatic potential at point r is equal to the electrostatic interaction energy between the unit positive charge placed at that point and the molecule 40,41 .According to Bader's suggestion, electrostatic potential is usually studied on the surface of van der Waals with an electron density of 0.001 a.u. 42,43Fig. 9 shows the surface diagrams of van der Waals colored by electrostatic potential of different systems and the surface penetration diagrams of van der Waals between anion and cation.Among them, different colors represent different electrostatic potential values, the redder the color, the more positive the electrostatic potential is, and the bluer the color, the more negative it is.The local minimum and maximum of electrostatic potential on van der Waals surface are drawn as cyan and orange spheres respectively.According to the principle of electrostatic potential complementarity, the positive region of electrostatic potential is easy to interact with the negative region 44 .It is not difficult to see from Fig. 9 that the positive potential region and negative potential region between m-DAPZ + and other anions penetrate each other, which reflects the complementary characteristics of electrostatic potential and also reflects the nature of electrostatic attraction interaction such as hydrogen bond and conjugated π-π.At the same time, this means that there is a strong attraction between anion and cation, which helps to form a stable composite system.The electrostatic potential on the upper surface of m-DAPZ + is distributed in a wide area, and the maximum and minimum values are about 141.16 kcal•mol −1 and 73.90 kcal•mol −1 respectively, which are the same in the three systems.In addition, the maximum values are basically near the amino group, which is mainly because the amino group is an electron-donating group, which makes the electrostatic potential on the surrounding surface appear positive.Correspondingly, the minimum values are all around the N atom on the pentazole, which is far away from the amino group and the electron density is relatively high.In the electrostatic potential diagrams of three anion systems, the minimum electrostatic potential is −5.34 kcal•mol −1 , −5.11 kcal•mol −1 and −5.82 kcal•mol −1 respectively, and the minimum electrostatic potential of N 3 − anion system is even smaller.In addition, the maximum point of its electrostatic potential is annular near the middle N atom, and its values are all −5.39 kcal•mol −1 .By comparing the van der Waals surface penetration diagram between ions with the whole van der Waals surface diagram, it can be found that the electron density of the whole ion system has changed greatly after the formation of the composite ion system.The values of extreme points in the positive potential region of m-DAPZ + van der Waals surface basically show a decreasing trend, which will mean that the electron density near m-DAPZ + increases.Compared with m-DAPZ + system alone, the polarity of some of its reaction sites decreases, and m-DAPZ + in the composite system will not easily react with external nucleophilic groups.However, the extreme points on the surfaces of the three anions van der Waals also show a decreasing trend, which shows that they are more likely to produce sites that react with external electrophilic groups.Generally speaking, m-DAPZ + is assembled with other anions, and they interact with each other and have strong attraction, which helps to stabilize the formation of the complex.

Energetic performance of assembled ion system
For a brand-new energetic compound with polynitrogen, its energetic performance is an eternal topic.The densities of these ionic compounds were calculated according to the methods proposed by Politzer and Rice (detailed information in SI) [45][46][47] .The enthalpy of formation 48,49  − ] is a new high energetic oxidizer which is very suitable to form a mixed system with negative oxygen balance explosives.What needs to be added is that in the calculation process, even though experts and scholars have supplemented some theoretical formulas in detail, in fact, these formulas still have some errors due to the limitations of the theory and the complexity of the composite system.Especially, when calculating the density, although Rice and others repeatedly consider the weak interaction between anion and cation, it is not difficult to see from the original text that the interaction of LP-π, π-π, etc. are not considered.It will mean that the volume calculated is still too large, which will lead to the decrease of the calculated density.On the contrary, its theoretical density will be greater, and of course, its energy will increase, the importance of its research is self-evident.Therefore, the appearance of the above three high energetic materials, especially [m-DAPZ + ][N 5 − ], will be very exciting.

Conclusion
In this study, DAPZ + and its assembled polynitrogen and high-energy ionic compounds were studied deeply and comprehensively by combining experimental means with theoretical calculation.Mass spectrometry analysis indicated that diamino-pentazole was first successfully prepared during the amination of [Na(H 2 O)(N 5 )] 2H 2 O by HOSA.The calculation of ground state energy and transition state shows that the reaction product is m-DAPZ + in solvent, which is a great breakthrough in the synthesis of pentazole derivatives.The outstanding performance of DAPZ + is perfectly reflected in the calculation of the assembled system.The results of AIM, IGMH analysis show that DAPZ + and N 5 − , N 3 − , C(NO 2 ) 3 − all have strong interactions, including hydrogen bonds, LP-π, π-π and extensive van der Waals interactions, which play a positive role in the formation of single crystals.In addition, SAPT results verify the above analysis and calculate the accurate interaction energy.As for the detonation performance of the three systems, the theoretical calculation results are exciting.The theoretical detonation velocities of the three systems are all above 9000 m•s −1 , especially [m-DAPZ + ][N 5 − ] reached 10,016 m•s −1 (EXPLO5 V6.05 value).Due to the limitations of theoretical methods, their actual detonation velocities will be higher.The specific impulse of [m-DAPZ + ][N 3 − ] can reach 319.8 s at a relatively low density.Even for [m-DAPZ + ][C(NO 2 ) 3 − ], it is a high-energy oxidizer that can be compared to HNF.These excellent properties may allow the above three kinds of polyanion systems to replace many energetic materials with poor performance at present and play a greater role in military, aerospace and other fields.Therefore, it has to be said that the synthesis of DAPZ + and the exploration of nitrogen-rich and high-energy compounds are relatively successful, which may promote the field of energetic materials to move towards polynitrogen.

Experimental preparation
All reagents and solvents were purchased from Aladdin and Energy Chemical as analytical grade and used without further purification.All experiments used deionized water with a resistivity of 18.2 MΩ.
For the cleavage experiments, these pentazoles reaction liquid were dissolved in CH 3 OH, and infused into the spectrometer's ion source at 10 μL•min −1 with a syringe pump.In tandem mass spectrometry (MS/MS) experiments, the desired negative-ion or positive-ion peaks are mass-selected and subjected to product ion mass analysis following collision-induced dissociation (CID) at variable collision voltages using N 2 as the collision gas.
Mass spectrometric ionization was carried out in electrospray ionization source (ESI) ionization mode with a source temperature of 250 °C.The detection was carried out on a TSQ Vantage triple quadrupole (QqQ) mass spectrometer (Thermo Fisher Scientific), equipped with an ESI turbo spray ionization source.The entire system was controlled via Xcalibur 2.2 software and data were processed using TraceFinder 3.1.

Figure 2 .
Figure 2. High-resolution mass spectra of the reaction solution.(a) Negative ion, mass spectrum of the reaction solution at 40 °C for 8 h.(b-d) Negative ion, full-range CID mass spectra of the mass selected, (m/z 100.945) peak recorded at collision voltages of −30, −120, and −200 V. (e) Positive ion, mass spectra of the reaction solution at 40 °C for 8 h.f positive ion, full-range CID mass spectrum of the mass selected, (m/z 102.128) peak recorded at collision voltages of 30 V.

Figure 6 .
Figure 6.Electron density topological diagrams of three ion systems (orange spheres represented bond critical points (BCPs).Bond paths were drawn as orange paths).

Figure 8 .
Figure 8. Decomposition diagram of interaction energy of three ions based on SAPT theory.

Table 2 . 5 +
Detonation performance of three systems and other energetic materials.a Calculated density.b Detonation velocity, values calculated using the K-J equation/EXPLO5 V6.05.c Detonation pressure, values calculated using the K-J equation/EXPLO5 V6.05.d Specific impulse of the neat compound using the EXPLO5 V6.05 program package at 7 Mbar chamber pressure.e Calculated heat of formation.f Oxygen balance (based on CO 2 ) for C a H b O c N d , 1600(c − 2a − b/2)/M W ; M W is the molecular weight.g Ref. 53 .h HNF: [N 2 H